Polymer size in dilute solutions in the good-solvent regime 



Sergio Caracciolo, Bortolo Matteo Mognetti 
Dipartimento di Fisica and INFN - Sezione di Milano I Universitd degli Studi di Milano 

Via Celoria 16, 1-20133 Milano, Italy 
e-mail: Sergio . CaraccioloSsns . it, Bortolo . MognettiOmi . inf n . it 

Andrea Pelissetto 
Dipartimento di Fisica and INFN - Sezione di Roma I 
Universitd degli Studi di Roma "La Sapienza" 

P.le A. Moro 2, 1-00185 Roma, Italy 
e-mail: Andrea . PelissettoSromal . inf n . it 

Abstract 

We determine the density expansion of the radius of gyration, of the hydro- 
dynamic radius, and of the end-to-end distance for a monodisperse polymer 
solution in good-solvent conditions. We consider the scaling limit (large de- 
gree of polymerization), including the leading scaling corrections. Using the 
expected large-concentration behavior, we extrapolate these low-density ex- 
pansions outside the dilute regime, obtaining a prediction for the radii for any 
concentration in the semidilute region. For the radius of gyration, compari- 
son with field-theoretical predictions shows that the relative error should be 
at most 5% in the limit of very large polymer concentrations. 

PACS: 61.25.Hq, 82.35.Lr 



I. INTRODUCTION 



For sufficiently high molecular weights, dilute and semidilute polymer solutions under 
good-solvent conditions exhibit a universal scaling behavior, as predicted by the renormal- 
ization group. For instance, if R is some typical dimension of the polymer, then R aN^, 
where N is the degree of polymerization, a is a constant that depends on density, temper- 
ature, and chemical details, and 1/ a universal exponent, u ~ 0.5876 (Ref. 6). From an 
experimental point of view, two quantities R are easily accessible in the dilute regime. The 
radius of gyration Rg is determined in small-angle scattering experiments. For small wave 
numbers q, the static form factor behaves as 

a/3 

that allows the determination of Rg. The second quantity of interest is the hydrodynamic 
radius Rh, which is related to the chain diffusion constant D that can be measured in 
inelastic scattering experiments:^ 

D = (1.2) 

where 77 is the viscosity. In theoretical work, another quantity is often used, the end-to-end 
distance Rg, which can, in principle, be accessible in experiments with polymers with tagged 
endpoints.^ 

In this paper we consider the density dependence of these three quantities. In general, if 
R is one of the quantities defined above or any functions of them, in the dilute regime we 
can write 

R{c) ^ R[l + Ric + R2C^ + 0{c^)], (1.3) 

where c is the polymer number density and R = R{c — 0). The coefficients Ri, R2, etc. 
depend on chemical details. However, os N ^ 00 the ratios Sn = Rn/Rf^ approach universal 
constants that are independent of temperature — as long as the polymer is in the good- 
solvent regime — and chemical properties. Thus, if N is large (scaling regime) we can write 
the universal expression 



R{c) 



= 1 + SlcRl 



+ s;{cRir + o[{cRin 



(1.4) 



R 



In many cases, polymers are not long enough to be in the scaling regime and thus it is 
important to include scaling corrections. As predicted by the renormalization group and 
extensively verified numerically, for N large but finite we have for the expansion coefficients 
of quantities related to the end-to-end distance and to the radius of gyration 



where A is a universal exponent whose best estimate is^ A = 0.515 ± 0.007iQ;QQo. The con- 
stants S* and Sn are universal. The constant k2 encodes all chemical details and nonuniversal 
properties of the solution — temperature, for instance. It can be determined by considering 
the scaling corrections to some renormalization-group invariant polymer quantity. We use 
the interpenetration radius ^ that can be determined from the osmotic pressure and rep- 
resents one of the most easily experimentally accessible quantities.^ In the large- limit it 
behaves as 



This relation defines the constant k2- The constant ^* is universal and known with very good 
precision:^" ^* = 0.24693±0. 00013. Other estimates are reported in Ref. 11. Universal ratios 
related to the hydrodynamic radius scale differently because of the presence of corrections 
proportional to N^~'^ (Ref. 12). Note finally that, given two different quantities Ri and R2 
that characterize the polymer size, their zero-density ratio R1/R2 approaches a universal 
constant RI2 for N 00. 

For the purpose of determining the universal constants defined above we perform a Monte 
Carlo simulation of the lattice Domb- Joyce model, considering walks of length varying be- 
tween 100 and 8000 and three different penalties for the self-intersections. Wc are able to 
obtain the leading density correction with a relative precision of 0.5% for Rg and Rg and of 
5% for Rh- We also estimate the second density correction, although with limited precision. 



Sfi — 




S: (1 + k2SnN-'^) , 



(1.5) 



* fa (1 + k2N-^) . 



(1.6) 
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These results show that density corrections are very small below the overlap concentration 
c* = 3/{4:7rRg), and that a simple approximation that assumes a linear dependence on the 
density works quite well up to c ~ c*. Once the virial expansion is known, we can try to 
resum it to obtain an interpolation formula that is valid in the whole semidilute regime. In 
Ref. 10 we have shown that a simple extrapolation of the virial expansion for the osmotic 
pressure that takes into account the expected large-concentration behavior works reasonably 
well. Comparison with experimental data and with other theoretical predictions shows that 
the error is less than 1% in the dilute region and rises at most to 5-10% deep in the semidilute 
regime. We use here the same strategy for the radii, determining simple expressions that 
reproduce the virial expansion for c — > and have the correct large-density behavior. Com- 
parison with other theoretical results shows that, at least for Rg, the interpolation formula 
has an error of at most 5%. 

The paper is organized as follows. In Sec. II we derive the virial expansion for a generic 
quantity that depends on the structure of a single polymer, like the radii defined above. 
In Sec. Ill we define the model and the quantites we compute. In Sec. IV we analyze the 
Monte Carlo results. In Sec. V we present our conclusions and compare our results with 
those obtained in other approaches. 



II. VIRIAL EXPANSION FOR SINGLE-POLYMER PROPERTIES 



We wish now to derive the density expansion of R{c), where i? is a single-polymer prop- 
erty. We consider L polymers in a volume V with configurational partition function 



Ql = 
d/i = 



exp 



(2.1) 
(2.2) 



where V^™*^^ is the sum of all terms of the Hamiltonian that correspond to interactions 
of monomers belonging to two different polymers i and j, ^A™*''^ is the contribution due 
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to interactions of monomers belonging to the same polymer i, and x^j is the position of 
monomer a belonging to polymer i. As usual, we define the Mayer function 

/. . = exp(-/?V-,f - 1. (2.3) 

Since i? is a single-polymer property, we can assume that it depends only on the coordinates 
Xq,i of the monomers of the first polymer. Thus, we wish to compute the expansion of 

(it:) = ^ y df, R{{^ai}) exp (^-P g - j , (2.4) 

for c — > 0. The term linear in the number density c = L/V is easily derived: 

l^^l.g) (/.)„ + ... (2.6) 

where 

Qo,L = j dii (e»)o = -^J diiO. (2.7) 
It follows that, for L, y — > oo at fixed concentration, we have 

{R)^R + c \ji{R) - RI^ , (2.8) 
where R is the zero-density value of R, and 

h = j d^ri2(/l2)0,ri„ (2.9) 
J,{R) = j dh^2 {Rf 12)0,r^,. (2.10) 

Here (^.r indicates an average over two independent polymers such that the first one starts 
at the origin and the second starts in r. In the mean values is a function of the coordinates 
of the first walk. 

The correction of order can be worked out analogously. We define 



5 



h = J dhi2dhn (/l2/l3/23)o,ri2,ri3 

J2{R)= J dh,2dh,, (i?/l2/l3)o,, 
MR) = J dh,2dh,^ (i?/l2/23)o,ri„r 
MR) = J d\,2dh,^ (i?/l2/l3/23)o,r.„r.3 ' 



',ri2,ri3 



,ri3 



(2.11) 
(2.12) 
(2.13) 
(2.14) 
(2.15) 



where (•)o,r2,r3 refers to averages over three polymers, the first one starting in the origin, the 
second in r2, etc. We obtain finally 



{R)^R + c Ji{R)-Rl2 
1 



+ 



r 



J2(R) + 2J3(R) + Ja{R) - SRTi - Rh - Ah^R) + Rlf\ + 0{^). (2.16) 



In the following we shall consider a lattice model for polymers. In this case, the previous 
expressions must be trivially modified, replacing each integral with the corresponding sum 
over all lattice points. 

III. MODEL AND OBSERVABLES 

Since we are interested in computing universal quantities, we can use any model that 
captures the basic polymer properties. For computational convenience we shall consider a 
lattice model. A polymer of length N is modelled by a random walk {fq, ri, . . . , r^} with 
\Ya — Ta+il = 1 On a cubic lattice. To each walk we associate a Boltzmann factor 



{)<a<P<N 



(3.1) 



with w > 0. The factor cr counts how many self-intersections are present in the walk. This 
model is similar to the standard self-avoiding walk (SAW) model in which polymers are mod- 
elled by random walks in which self-intersections are forbidden. The SAW model is obtained 
for w — -\-oo. For finite positive w self-intersections are possible although energetically pe- 
nalized. For any positive this model — hereafter we will refer to it as Domb- Joyce (DJ) 
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model — has the same scaUng Hmit of the SAW model^'^ and thus allows us to compute the 
universal scaling functions that are relevant for polymer solutions. The DJ model has been 
extensively studied numerically in Ref. 8. There, it was also shown that there is a particular 
value of w, w*, for which corrections to scaling with exponent A vanish: the nonuniversal 
constant k2 is zero for this particular value of w. Thus, simulations at w — w* are partic- 
ularly convenient since the scaling limit can be observed for smaller values oi N. In Ref. 8 
w* was shown to be approximately equal to 0.505838 {e~^* — 0.603), while in Ref. 10 it was 
found w* ^ 0.47 ±0.02. 

The DJ model can be efficiently simulated by using the pivot algorithm. For the 
SAW an efficient implementation is discussed in Ref. 18. The extension to the DJ model is 
straightforward, the changes in energy being taken into account by means of Metropolis test. 
Such a step should be included carefully in order not to loose the good scahng behavior of 
the CPU time for attempted move. We use here the implementation discussed in Ref. 19. 

We shall be interested in the quantities that characterize the size of the polymer and 
have been defined in the introduction. In the lattice model they are defined as follows: 



p2 _ 

^ 2{N 



(3.) 



+ 1)^ |r. - 

Rl = {(r„ - r„)^> . (3.4) 

The expansion coefficients Sn have been determined by using the expressions reported in 
Sec. II with 

f.j = e-^^'i - 1 = 6-"""'^ - 1 (3.5) 

a/3 

Here r^j is the position of monomer a of polymer i. To compute the corresponding lattice 
sums, we use a simple generalization of the hit-or-miss algorithm discussed in Ref. 20 and 
in the appendix of Ref. 10. 



In the dilute limit we write 

= 1 + S,,,{cRl) + S,,,{cRir + ■■■ (3.7) 

f^l + S,,,{cRl) + S2,e{cRy + --- (3.8) 

Re 

The coefficients -S'^^:;^ have a finite hmit for — > oo. Including the leading scaling correction, 
in the scahng hmit they behave as follows: 

^n,#(A^) = 5;#(l + ^ + •••), (3.9) 

where the constants S*^^ and are universal. The factor k2 depends instead on chemical 
details, temperature, etc. and is fixed by the large- A?" behavior (1.6) of the interpenetration 

radius 

Beside the coefficients parametrizing the density dependence of the radii, we can also 
consider the large- A?" behavior of the radii that we parametrize as 

Rl = a,Ar2- ^1 ^ A;2r,Ar-A + (3.10) 

Rl = GeN^" (1 + k2reN^^ + •■■). (3.11) 

The constants and are universal as well as the ratio A*^ = ag/ue- This quantity can 
also be determined directly by considering 

= |f ^ A;^ (1 + /c^r^eAT-^ + ...), (3.12) 

where r^e ^ Tg - r^. 

Quantities involving the hydrodynamic radius must be analyzed differently, due to the 
appearance of new nonuniversal corrections proportional to N"~^. For instance, as discussed 
in Ref. 12, 1/Rh, which is the quantity measured in MC calculations, scales as 

^^'^X, + a,N--) + l^, + - = Z(l + ^N^-^+a,N-- + :\ . (3.13) 



Rh N N" \ ai 

The analytic term gives rise to a new correction proportional to A^''"^ — _/V~°'^^^. It de- 
cays slower than the standard one proportional to A^""^ = 7V~°-^ and does not satisfy any 
universality property. We also consider the density expansion 

8 



9^ = 1 + S^,H{cRl) + S,,„{cRlf + • • • (3.14) 
tin 

and the ratio 

A,H ^ 1^. (3.15) 
Rh 

For iV — > oo, Sn,H and converge to universal constant 5* and A*jj. Note that it is 
more convenient to consider Rg/Rn than Rn/Rg, since the latter quantity has additional 
corrections proportional to N'^'^"^ — 7V~°-^^^, etc. 



IV. ANALYSIS OF THE MONTE CARLO DATA 

We have performed three sets of simulations ai w — 0.375, 0.505838, 0.775, using walks 
with 100 < < 8000. Results are reported in Tables I, II, III, and IV. The numerical 
data that concern R^ and R^ have been analyzed as discussed in Refs. 21, 10. Therefore, if 
Q — Age, Si^g, S2,g, Si^e, and S2,e, we fit the data to 

Q(iV,.).«.^^^£^. (4.1) 

For A we use the best available estimate: A = 0.515±0.007lo:oo2 (Ref. 8). The term l/AT^^.es 
should take into account analytic corrections behaving as 1/N, and nonanalytic ones of the 
form N'"^^, 7V~^2 (^/\^ jg ^j^g next-to- leading correction-to-scahng exponent). As discussed 
in Ref. 21, one can lump all these terms into a single one with exponent A2,eff = 1.0 ± 0.1. 
The results are reported in Table V. The reported error is the sum of the statistical error and 
of the systematic one due to the uncertainty on A and A2,eff- In order to detect additional 
corrections to scaling we have repeated the fit several times, including each time only the data 
with N > A^min- More precise estimates can be obtained by using the universality of the ratios 
of the correction-to-scaling amplitudes. As we did in Ref. 10, we can simulteneously analyze 
two quantities Qi and Q2, using the fact that aQ-^{w) = ri2a(32('"^)) with ri2 independent of 
w. Thus, one can fit simultaneously all data for Qi and Q2, taking as free parameters Ql, 
Q2, aQ^{wi), CQ^{wi), CQ^{wi), and ri2- We have reanalyzed the data for Si^g, S2^g, S'l^g, and 



5*2,6 using this method, taking always as second variable Rg/Rl- The results arc reported in 
Table VI. The improvement is quite significant, the errors decreasing by a factor of 2-3. We 
do not observe significant residual scaling corrections and take the result with A^min = 500 
as our final estimate. Therefore 

= 0.15988 ± 0.00004, (4.2) 
= -0.3152 ± 0.0007, (4.3) 
^2*^ = -0.087 ± 0.005, (4.4) 
= -0.3853 ± 0.0007, (4.5) 
= -0.027 ± 0.006. (4.6) 

Note that all density corrections are negative — this is not surprising since, by increasing 
the density, the polymer size decreases. Moreover, both for Rg and Rl, the second density 
correction is significantly smaller than the first one. There are in the literature several 
estimates of A*^ (see Ref. 11 for an extensive list of references). The most precise ones are: 
A*g^ ^ 0.1599±0.0002 (Ref. 20) and A^^ = 0.15995±0.00010 (Ref. 22). They are both in good 
agreement with our result. We can also compare the estimates of S*g and 5'2 ^ with those 
obtained by using one-loop perturbation theory:^'^^ 5"*^ ~ —1.19 and 5*2 ^ ~ 5.74. They 
both significantly differ from our estimates, indicating that the perturbative renormalization 
group is not accurate (at one-loop order) for these quantities. 

In the analysis of quantities involving the hydrodynamic radius we should take into ac- 
count the corrections proportional to A^'^"^ — iV~°-^^^. These new corrections are particularly 
strong and depend only slightly onw, as it can be seen from the results in Tables I, IV: for all 
values of w, Rg/Rn and Sn,H show the same monotonic behavior, indicating that corrections 
proportional to are much smaller than those proportional to A^''"^ (rememeber that 

the corrections approximately vanish for w — 0.505838 and have opposite sign in the 
other two cases). In order to determine the scaling-limit value of AgH, Si^h, and S2,h we fit 
the data to 
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«iV.») = «- + ^ + ^. (4.7) 

Results are reported in Table V for several values of A^min- As before, the error takes into 
account the uncertainty on u and A. In all cases, the results show small trends with A^min 
indicating that there are some small additional scaling corrections that are not taken into 
account by the fit Ansatz (4.7). We quote as our final results: 

= 1.5810 ±0.0010 (4.8) 
SI jj = 0.082 ± 0.004 (4.9) 
Sl^H = 0-050 ± 0.020. (4.10) 

The error should be large enough to take into account all residual systematic corrections. 
Note that also in this case density corrections are small. The result for A*jj is in good 
agreement with that obtained in Ref. 12: A*^ = 1.591 ± 0.007. 

Finally, we compute the universal scaling-correction coefficients: our data are precise 
enough to provide estimates of Si^g, Si^^i ^ge; ^e; and Tg. The first three quantities can be 
determined as discussed in Ref. 21. If Q{N, w) is a renormalization-group invariant quantity 
we consider 



\l>* 



Q{N,w,)-Q{N,W2) 



(4.11) 



which should converge asymptotically to the correction-to-scaling amplitude with corrections 

of order N'^"^ , with Aeff = 0.5 ± 0.1. In practice we use the data with wi — 0.375 and 
W2 = 0.775 and fit Tq{N) to 

We find that the correction-to-scaling coefficients are quite small: 

Vge = -0.0040 ± 0.0004, (4.13) 

s^ g = -0.050 ± 0.010, (4.14) 

= -0.050 ±0.020. (4.15) 
11 



As for S2,g and S2,e we only obtain upper bounds: S2,g = —0.1 ± 0.2, |s2,e| ^ 0.2. In order to 
determine Tg and Vg we first define Qr{N) = R'^{2N)/R'^{N) and then 



Tr{N) = 



Qr{N,Wi)-Qr{N,W2] 



(4.16) 



1-2^ [^{2N,wi) -^{2N,W2), 
that converges to Tg or r^, depending on the radius one is considering, with corrections of 



order N '^<'^. We obtain: 



rg = -0.36 ± 0.07, (4.17) 
re = -0.28 ± 0.08. (4.18) 

Finally, if we analyze as before the quantity 

T (m - Qn,{N,wi)-QR^{N,W2) 

QrAN,Wi) -Qr^{N,W2) 

we obtain a direct estimate of rg/ve'- 

r-g/r-e = 1.2 ± 0.3. (4.20) 

This is consistent with the estimate 

!l = 1 + !^ = 1.014 ±0.004, (4.21) 
re re 

obtained by using r^g and r^. These results can be compared with those appearing in the 
hterature. Ref. 24 gives —0.5 < < —0.4 and Vg/ve — 1.25 ± 0.04, while Ref. 8 gives 
Ve — —0.41 ± 0.01. The estimates of Vg and are consistent with ours, while that for r^/rg 
differs significantly, given the small error bars. 



V. CONCLUSIONS 

The results of the previous Section allow us to determine the radii in the dilute regime. 
We find 
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-J-^l- 0.0752$p - 0.0050$J + • • • + A;$(0.003$j, + hg^l ■■■), (5.1) 

^ = 1 - 0.0920$p - 0.0015$J + • • • + A;$(0.004$j, + h^^l • • •), (5.2) 

^ = 1 + 0.0195$p + 0.003$^ + • • • , (5.3) 
Rh 

where we have introduced the polymer packing fraction 

^ 3 3 ^ ' 

N A is the Avogadro number, M the molar mass of the polymer, c and p the number density 
and the ponderal concentration respectively. The constants hg and 6e depend on si^g and S2,e 
and thus we only have an upper bound: \bg\ < 1 • 10~^, |6e| ^ 2 • 10""^. The constant A;$ takes 
into account scaling corrections. It scales as N'^^ where N is the degree of polymerization, 
and therefore vanishes in the scaling limit. It is defined by the small-density behavior of the 
osmotic pressure:-'^^ 

Z = ^ « 1 + 1.313$p + 0.559$J - 0.122$J + . . . + k^{^^ + 1.13$J + •••)• (5.5) 

Expansions (5.1), (5.2), (5.3) are supposed to be valid up to $p ~ 1. In this range, the linear 
term in the density dominates, the next-to-leading term being quite small, at variance with 
the behavior of the osmotic pressure. Note also that the density dependence is very small in 
all cases. The nonuniversal parameter allows us to take into account the deviations from 
the scaling limit due to the finite degree of polymerization (as long as scaling corrections are 
small). It can be obtained from measurements of the osmotic pressure by using Eq. (5.5). 
Once is known, everything is fixed and the previous expressions give us the radii in the 
whole dilute regime $p < 1. Note that corrections to scaling in the radii are much smaller 
than in the osmotic pressure. 

It is wortwhile to try to extrapolate the above-reported expansions outside the dilute 
regime. As in Ref. 10 we will use a simple interpolation formula that takes into account the 
large- $p behavior of the radii. For $p — > oo we have^ 
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5^ ^ $-(2v-l)/(3i.-l) ^ $-0,230 gN 
^2 P ' ^ ■' 

which is a consequence of the fact that, at fixed nonvanishing monomer concentration, poly- 
mers behave as Gaussian coils as N ^ oo. It is easy to write down interpolation formulas 
that take into account Eq. (5.6) and reproduce the previous expansions for — > oo: 

^ = (1 + 0.655*„ + 0.28*^)-°•^^^ 

d2 

^ = (1 + 0.801*. + 0.37*^)-°•^^^ 



R 



^ = (1 + 0.34$j, + 0.10*2)0 0574_ ^5 

We do not attempt to include also scaling corrections here, since we have only estimated 
one single coefficient. Of course, these expressions are very precise in the dilute limit. In 
order to assess their precision in the semidilute regime we can compare with field-theoretical 
estimates^^"^^'^ and with Monte Carlo results. In Fig. 1 we compare expression (5.7) for 
the radius of gyration with the expression reported in Ref. 25. The two curves are quite 
close. They differ at most by 3.2% (for $p ^ 1.5) and, as $p oo, they give R^/Rl ~ 
1.18$;°-230 (Ref. 25) and R^/Rl ^ 1.16$;°-230 [Eq. (5.7)]. Note that the largest differences 
are observed for small values of $p (for $p = 0.4 the curves differ by 2%) where our expression 
is by construction very accurate. This is related to the fact that, as observed in Sec. IV, 
the low-density coefficients for the radii are poorly determined by one-loop field-theoretical 
calculations. Similar results are obtained by using the expressions reported in Ref. 5. For 
instance, they predict R^/Rl ~ 1.14$-°-230 [see Eq. (18.12) in Ref. 5]. In Fig. 1 we also 
report some Monte Carlo data taken from Ref. 28. They are somewhat lower (differences 
are less than 5%) than prediction (5.7) and in better agreeement with the field-theoretical 
curve, but they could be affected by finite-size and scaling corrections. This comparison 
shows that our expression for the density dependence of Rg is reasonably accurate, the error 
being apparently less than 5% in the whole semidilute region. We are not aware of similar 
field-theoretical results for the other radii, so that we are unable to estimate the precision of 
the extrapolation for Re and Rh- 
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TABLES 



TABLE I. Estimates of the ratios Age = Rl/Rj and AgH = Rg/Rn- 



N 


w = 0.375 


w = 0.505838 


w = 0.775 


Age 


100 


0.161771(15) 


0.1609797(95) 


0.159965(14) 


250 


0.160842(15) 


0.160283(11) 


0.159637(10) 


500 


0.160487(15) 


0.160079(11) 


0.159612(14) 


1000 


0.160293(21) 


0.160017(16) 


0.159644(24) 


2000 


0.160162(22) 


0.159910(15) 


0.159691(27) 


4000 


0.160067(22) 


0.159917(17) 


0.159753(23) 


8000 


0.160060(27) 


0.159888(19) 


0.159759(26) 




100 


1.246920(42) 


1.253683(42) 


1.263604(39) 


250 


1.348559(45) 


1.354758(52) 


1.363108(43) 


500 


1.404691(68) 


1.410316(71) 


1.417575(62) 


1000 


1.44766(10) 


1.452355(91) 


1.458507(99) 


2000 


1.48008(15) 


1.48402(15) 


1.48946(16) 


4000 


1.50436(23) 


1.50827(22) 


1.51251(20) 


8000 


1.52336(30) 


1.52607(32) 


1.52933(32) 
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TABLE II. Estimates of the coefficients Si^g and 6*2,3 radius of gyration. 



N 


w = 0.375 


w = 0.505838 


w = 0.775 




100 


—0 28885('58'l 


—0 30737('43'l 


—0 32464('55'l 






U.O±00-Ll J 




500 


-0.30663(71) 


-0.31391(45) 


-0.32202(65) 


1000 


-0.3093(10) 


-0.31533(72) 


-0.3184(11) 


2000 


-0.3121(11) 


-0.31562(74) 


-0.3183(10) 


4000 


-0.3105(11) 


-0.31412(77) 


-0.3161(13) 


8000 


-0.3144(13) 


-0.31629(87) 


-0.3166(14) 




100 


-0.0451(49) 


-0.0783(45) 


-0.1324(55) 


250 


-0.0646(59) 


-0.0905(48) 


-0.1234(47) 


500 


-0.0752(65) 


-0.0787(51) 


-0.1190(72) 


1000 


-0.071(10) 


-0.0950(73) 


-0.082(11) 


2000 


-0.059(11) 


-0.0920(78) 


-0.105(10) 


4000 


-0.079(11) 


-0.0883(81) 


-0.080(13) 


8000 


-0.099(11) 


-0.0896(89) 


-0.080(12) 
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TABLE III. Estimates of the coefficients Si^e and S2,e for the end-to-end distance. 



N 


w = 0.375 


w = 0.505838 


w = 0.775 


Sl,e 


100 


—0 ^M14(79) 


—0 37943('55'l 


— n 4n43Q('68'l 




U.OU 1 oo\ i Oj 






500 


-0.37364(85) 


-0.38442(56) 


-0.39570(91) 


1000 


-0.3780(13) 


-0.38556(85) 


-0.3911(11) 


2000 


-0.3814(13) 


-0.38608(91) 


-0.3891(13) 


4000 


-0.3801(13) 


-0.38397(92) 


-0.3866(14) 


8000 


-0.3848(16) 


-0.38609(97) 


-0.3878(16) 


S2,e 


100 


0.0050(67) 


-0.0271(62) 


-0.0862(86) 


250 


-0.0058(86) 


-0.0364(59) 


-0.0694(65) 


500 


-0.0139(82) 


-0.0225(69) 


-0.062(10) 


1000 


-0.009(13) 


-0.040(10) 


-0.025(13) 


2000 


-0.006(15) 


-0.042(10) 


-0.038(14) 


4000 


-0.040(14) 


-0.018(10) 


-0.008(14) 


8000 


-0.039(16) 


-0.035(11) 


-0.025(14) 
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TABLE IV. Estimates of the coefficients Si^h and S2,h for the hydro dynamic radius. 



N 


w = 0.375 


w = 0.505838 


w = 0.775 




100 


0929^(19) 


09747('18'l 


10014('19') 




U.Uc/ZpOc/I ZpU J 


nQ4Qi (9^ \ 


U.Uc/O 1 0\^\J J 


500 


0.09121(25) 


0.09221(28) 


0.09313(20) 


1000 


0.08991(39) 


0.09031(40) 


0.09097(34) 


2000 


0.08739(48) 


0.08874(53) 


0.08768(50) 


4000 


0.08708(68) 


0.08822(53) 


0.08747(59) 


8000 


0.08631(84) 


0.08612(73) 


0.08564(69) 


S2,H 


100 


0.0668(19) 


0.0851(21) 


0.0999(24) 


250 


0.0731(23) 


0.0826(22) 


0.0967(23) 


500 


0.0782(30) 


0.0822(27) 


0.0902(29) 


1000 


0.0811(42) 


0.0792(40) 


0.0929(41) 


2000 


0.0730(50) 


0.0911(62) 


0.0747(62) 


4000 


0.0842(78) 


0.0819(72) 


0.0748(81) 


8000 


0.0773(97) 


0.0738(94) 


0.058(10) 
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TABLE V. Fit results for several values of iVn 



Nrr 



100 



Nrr 



250 



Nrr 



500 



Nrr 



1000 



A 



ge 



A 



gH 



S. 



2,3 



'-'1,6 



"2,6 



S; 



2M 



0.15991(4) 
1.5831(4) 

-0.3144(8) 
-0.082(6) 

-0.3846(9) 
-0.023(8) 
0.0809(7) 
0.068(7) 



0.15989(4) 
1.5819(4) 
-0.3150(10) 
-0.090(9) 
-0.3851(11) 
-0.031(12) 
0.0815(10) 
0.059(12) 



0.15988(4) 
1.5811(7) 
-0.3152(15) 
-0.093(15) 
-0.3851(18) 
-0.032(19) 
0.0821(18) 
0.044(22) 



0.15995(7) 
1.5806(13) 
-0.3165(30) 
-0.095(27) 
-0.3873(36) 
-0.036(36) 
0.0866(36) 
0.050(44) 



TABLE VL Eslmiatcw of the density corrections lo Ry and i?, . Combiued fits with data for Ay,.. 

iVn,in = 100 iV^i^ = 250 Ar^i^ = 500 AT^in = 1000 

Slg -0.3149(9) -0.3151(6) -0.3152(7) -0.3143(13) 

Slg -0.087(5) -0.086(4) -0.087(5) -0.087(9) 

-0.3846(13) -0.3852(8) -0.3853(7) -0.3843(14) 

-0.028(8) -0.027(5) -0.027(6) -0.034(15) 
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FIG. 1. Plot of Rg/Rg versus <I>p. We report expression (5.7), "V", the field-theoretical predic- 
tion of Ref. 25, "FT", and the Monte Carlo data of Ref. 28. 
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